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Investigation  of  Aeroelastic  Flow  Control  of  a  Fluttering  Wing 
with  HPCMP  CREATE™-AV  Kestrel 


Casey  Faglcy* *,  Jurgen  Seidel*,  Thomas  McLaughlin* 
Department  of  Aeronautics,  U.S.  Air  Force  Academy,  CO  80840,  USA 


The  aeroelastic  behavior  of  a  finite  aspect  ratio  (At-  6)  NACAO018  wing  is  computationally  an- 
alyzed.  HPCMP  CREATE1  M-AV  Kestrel,  a  fully  coupled  computational  fluid  dynamic  (CFD)  and 
computational  structural  dynamic  fCSD)  code,  is  used  to  compare  tile  effect  of  blowing  on  a  rigid 
and  an  aeroelastic  ally  deforming  wing.  Externally  controlled  blowing  slots  distributed  along  the  span 
of  the  wing  are  used  to  inject  mass  into  the  flow  field  to  achieve  open-loop  flow  control.  The  results 
indicate  that  control  by  blowing  has  a  pronounced  effect  on  the  flow  field  and  therefore  on  the  aerody¬ 
namic  coefficients.  For  the  rigid  wing,  the  lift  is  increased,  as  are  the  pitching  and  rolling  moments. 
When  aeroelastic  deformation  is  considered,  the  picture  changes  significantly.  First,  the  flexible  wing 
showed  higher  lift  and  drag  compared  to  the  rigid  wing  due  to  the  increased  local  angle  of  attack  in 
the  outboard  section  of  the  wing.  Furthermore,  the  pitching  and  rolling  moments  were  significantly 
reduced. 
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1.  Introduction 

Af rq e l a st ic it y  ls  a  very  well  studied  scientific  field  combining  structural  dynamics,  aerodynamics  and  inter- 
actions  thereof.  Typically,  aeroelasticity  is  represented  by  Collars  triangle  of  forces  which  shows  the  interplay 
between  forces  present  in  an  aeroelastic  system.  A  simplified  triangle  offerees  is  depicted  in  Figure  I , 

The  transfer  of  energy  between  inertial,  clastic,  and  aerodynamic  forces  produces  a  potentially  unsteady  aero¬ 
elastic  behavior.  This  unsteady  fluid-structure  interaction  is  a  prominent  design  consideration  for  current  air  vehicles: 
ior  instance,  at  certain  operating  conditions  (i.e.  Mach  number,  altitude,  fuel  levels,  cargo  weight,  etc.)  aeroelastic 
instabilities  may  be  present  and  be  detrimental  to  the  aircraft  structure.  This  instability  is  ty  pically  identified  as  flut¬ 
ter,  or,  in  the  situation  of  a  marginally  stable  system,  a  limit  cycle  oscillation.  Aircraft  design  needs  to  maintain  a 
reasonable  safety  margin  from  any  aeroelastic  instability;  therefore,  a  wide  variety  of  methods  exist  to  predict  flutter 
margins  [  -  ].  Each  of  these  methods  have  one  unified  goal  which  is  to  identify  the  true  damping  of  the  aero-structure 
system.  This  damping  gives  the  stability  criterion  as  a  function  of  operating  condition  (mainly  Mach  number  and 
altitude).  Current  aircraft  designs  are  pushing  these  stability  margins  to  achieve  lighter  weight  and  therefore  more 
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efficient  flight  vehicles;  thus,  the  need  to  understand  the  source  of  these  aeroelastic  instabilities  and  potentially  be  able 
to  suppress  the  instability  on  order  to  expand  the  flight  envelope  is  greater  than  ever. 

The  typical  approach  to  suppressing  aeroelastic  instabilities  is  by 
using  control  surfaces,  be.  flaps  or  ailerons,  to  counter  act  unsteady 
aerodynamic  loads  due  to  structural  displacements.  In  terms  of  ac¬ 
tuation  most,  if  not  all,  active  aerodasttc  control  endeavors  utilize 
a  conventional  trailing  edge  flap  to  alter  the  aerodynamics.  For  in¬ 
stance,  in  the  case  of  the  design  of  the  X-56 A,  a  NASA  designed 
Multi-Utility  Technology  Testbed  for  flutter  suppression  techniques, 
the  trailing  edge  flap  actuator  bandwidth  was  10  Hz  [  j  and  the  flut¬ 
ter  instabilities  were  on  the  order  of  5-15  ilz  [  ]t  Alternatively,  sen¬ 
sors  such  as  accelerometers  and  strain  gauges  have  long  been  the 
conventional  devices  of  choice  to  instantaneously  measure  the  wing 
deformation.  Both  sensors  rely  upon  the  structural  deformation  for 
sensing  of  the  aeroelastic  phenomena.  The  bandwidth  demand  for 
flutter  suppression  is  a  huge  concern  for  the  stability  and  robustness  of  the  control  system  and,  ultimately,  leads  to 
the  motivation  of  this  work:  to  employ  active  flow  control  to  augment  the  local  flow  features  to  suppress  aeroelastic 
instabilities. 

Active  flow  control,  which  introduces  a  small  amount  of  energy  to  the  flow  (e,g,  by  blowing  fluid  into  the  main 
flow),  allows  for  direct,  local,  and  instantaneous  modification  of  the  flow  field.  The  perturbation  of  the  flow  is  amplified 
through  a  natural  fluidic  instability  which  affects  the  global  flow  state.  Thus,  the  control  effort  remains  small  by 
exploiting  the  fluidic  instability  and  the  bandwidth  of  the  system  is  equivalent  to  the  convective  frequency  of  the 
flow  [  ].  In  the  past,  both  passive  and  active  (including  open-  and  closed-loop)  flow  control  techniques  have  been 
exhaustively  studied  (see  review  articles,  e.g.  [  ,  ]),  but  upon  a  survey  of  the  literature,  implementation  of  active 
closed-loop  flow  control  to  real-life  applications  has  largely  eluded  the  research  and  engineering  communities  to 
date  [  ],  The  development  of  a  functioning  control  system  for  a  range  of  flow  conditions  has  proven  to  be  a  formidable 
challenge. 

The  research  presented  in  this  paper  utilizes  Kestrel,  a  fully  coupled  computational  fluid  dynamic  (CFD)  and  com¬ 
putational  structural  dynamic  (CSD)  code,  to  assess  the  aeroelastic  behavior  of  a  generic  wing.  Externally  controlled 
mass  injection  slots  are  used  to  augment  local  flow  features,  and  the  global  response  is  characterized.  It  is  shown 
through  these  high  fidelity  simulations  that  flow  control  on  aeroelastically  deforming  surfaces  is  possible.  Also,  a  flow- 
field  analysis  is  performed  and  interesting  relationships  betw  een  forcing  and  structural  responses  are  obtained, 

II.  Setup 

Simulations  w  ere  performed  on  a  cantilevered,  rectangu¬ 
lar,  flexible  NACAG0I8  wing  with  a  span  of  b  —  0,6m, 
a  chord  of  c  —  0.1m,  and  an  aspect  ratio  of  The 

geometry  is  shown  in  Fig.  2;  note  that  the  origin  of  the 
coordinate  system  is  at  the  root  leading  edge,  with  the 
x-axis  pointing  downstream,  the  y-axis  in  the  spanwise 
direction,  and  the  z-axis  upward. 

Four  rectangular  blowing  slots  were  integrated  on  the 
upper  surface  near  the  leading  edge  of  the  wing.  The 
slots  were  placed  in  1/2  chord  increments  from  the  tip 
toward  the  root  of  the  wing.  The  slot  dimensions  are 
l  mm  by  10mm.  The  orientation  of  the  slots  was  chosen 
to  be  angled  upward  by  30E?  with  the  slot  blowing  surface 
on  the  chord  line  at  the  10%  chord.  The  slot  location 
and  geometry  are  also  shown  in  Fig,  2{bottom).  Since 
the  slots  were  placed  in  an  area  of  a  favorable  pressure 
gradient,  it  is  expected  that  the  effect  forcing  the  flow 
by  blowing  out  through  the  slots  will  be  most  effective 
at  large  angles  of  twist  (i,e„  large  local  angles  of  attack) 
where  the  flow  is  dose  to  separation  near  the  leading 

Fi  g u  rc  2,  N  ACA00 1 8  geometry  a  ltd  blowing  slot  co  ttfigu  r  a  tfon. 

The  free  stream  velocity  was  set  to  £/«  =  34m/s,  cor- 


- .  . - . -  o.io— 

•  30,00" 


Figure  1.  Modified  Collar's  triangle  representing  the 
forces  at  play  In  an  aeroelastic  system. 
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responding  to  a  Mach  number  of  M  =  0.  J .  The  chord  based  Reynolds  number  is  Re^  =  1 60, 000  and  the  angle  of  attack 
was  set  to  a  —  6  deg. 

A*  Coupled  Aeroelastlc  CFD 

The  HPCMP  CREATE™-AV  Kestrel  software  suite  was  employed  for  the  simulations  in  this  research.  The  software 
suite  is  designed  to  allow  for  development  of  extensions  such  as  a  flow  control  component,  [  -  ]  The  underlying 

program  structure  is  based  upon  an  event  based  model,  where  the  various  components  in  the  simulation  subscribe  to 
and  trigger  events*  This  architecture  allows  for  great  flexibility  of  simulations  without  sacrificing  performance. 

Kestrel  is  based  on  the  AVUS  flow  solver  developed  at  the  Air  Force  Research  Laboratory  at  Wright- Patterson 
Air  Force  Base.  It  is  a  multidisciplinary,  second-order  in  space  and  time,  cell -centered,  finite-volume  Navier- Stokes 
solver  for  grids  with  arbitrary  cell  topologies.  Kestrel  features  different  flux  schemes,  limiters  and  turbulence  models. 
Aeroelastic  analysis,  moving  control  surfaces,  overset  grids  and  rigid-body  motions  are  some  of  the  features  in  the 
current  version  of  Kestrel. 

Kestrel  v5. 1  features  a  number  of  turbulence  models,  including  Spalart-A!lmaras(SA)  [  ■],  Spalart-AHmaras  with 
rotation  correction  (SARC)  [  ],  Mentor's  model  (Menter)  [  ],  and  Mentor's  shear  stress  transport  (SST)  \  ];  all 

RANS  models  can  be  operated  in  Delayed  Detached-Eddy  Simulation  (DDES)  [i  j  mode  to  provide  hybrid  RAN- 
S/LES  (Large  Eddy  Simulation)  simulation  capability.  The  simulations  presented  in  this  paper  were  performed  using 
the  DDES  approach  with  the  SA  turbulence  model. 

To  accomplish  the  final  research  goal  of  this  project,  namely  feedback  flow  control  of  a  flexible  wing,  the  capability 
of  implementing  externally  controlled  boundary  conditions  using  the  Kestrel  Software  Development  Kit  (SDK)  is 
crucially  important.  Feedback  flow  control  is  accomplished  in  general  via  sensors  in  the  flow,  a  controller,  and  an 
actuator.  A  generic  'ComponentX’  is  provided  with  the  SDK  and  was  modified  to  implement  feedback  flow  control, 
hi  the  current  implementation,  sensor  data  are  provided  using  the  access  to  grid  and  solution  quantities  provided 
within  the  SDK  through  the  'Data  Warehouse’  infrastructure.  [  ]  The  controller  can  be  implemented  directly  in 

C  omponentX' ;  the  output  of  the  controller  is  then  written  back  to  the  respective  mass  flow  source  boundary  condition. 
Wide  this  procedure  to  set  the  boundary'  condition  was  used  for  the  current  simulations,  results  shown  tn  this  paper 
do  not  yet  utilize  the  full  feedback  mechanism;  instead,  the  mass  source  boundary  condition  is  set  to  a  fixed  values. 
Therefore,  the  forcing  is  referred  to  as  ’open-loop'  control  rather  than  'closed-loop’  or  feedback  control.  While  this 
type  of  control  could  also  be  achieved  using  the  standard  mass  flow  source  boundary'  condition,  it  was  important  to 
implement  it  in  this  manner  to  test  the  infrastructure  for  eventual  implementation  of  feedback  control* 

The  geometry  in  Fig*  J  was  created  in  Sotidworks®  and  imported  into  Simcenter  Solid  Mesh*  [  ]  The  geometry 
was  completed  by  a  symmetry  plane  at  the  wing  root  and  a  hemispherical  farfield  with  a  radius  of  10  times  the  wing 
span.  T  he  grid  spacing  on  the  wing  w  as  chosen  to  match  the  grids  which  were  developed  in  the  companion  paper  j_  J, 
which  includes  a  spatial  and  temporal  resolution  study.  The  wall  normal  spacing  resulted  in  y+  <  1  for  the  whole  wing 
expect  inside  the  blowing  ports,  where  yH  <  3  for  the  case  w  ith  blowing;  the  complete  grid  had  20*3  million  cells.  The 
w  ing  geometry  with  the  grid  on  the  symmetry  plane,  the  w  ing  surface  grid  and  a  close-up  of  the  grid  at  a  blowing  slot 
are  shown  in  Fig.  3. 


a)  N  AC  ADO  IS  Wing  t>)  Surface  Mesh  c)  Blowing  Slot 


Figure  3*  20.3  million  cell  grid  of  the  Unite  aspect  ratio  NACA0018  wing* 


R,  Structural  Model 

To  utilize  the  fluid-structure  interaction  components  within  Kestrel,  a  modal  structural  model  of  the  geometry  is  neces¬ 
sary.  For  the  basic  problem  of  a  linearly  elastic  wing,  where  the  material  obeys  Hooke’s  Law,  the  system  of  equations 
to  be  analyzed  is  a  second  order  ordinary  differential  equation  describing  a  mass-spring-damper  system*  This  system 
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can  be  written  after  a  Laplace  transform,  as 

[Ms2  +C$+K]q  =  F,  (]) 

where  M  is  the  mass  matrix,  C  is  the  damping  matrix,  K  is  the  stiffness  matrix,  q  is  the  vector  of  degrees  of  freedom 
(6-DoF  for  each  node,  i,  three  translations  and  three  rotations),  5  is  the  Laplace  variable,  and  F  is  the  force  acting  at  the 
nodal  location.  For  practical  implementation,  the  geometry  is  discretized  and  a  finite  element  approach;  the  resulting 
mass,  damping  and  stiffness  matrices  are  formed  and  the  system  of  equations  are  formulated  as. 


A 

"  0 

1  ' 

[ql 

A 

K 

.  M 

-c 

An  eigen-analysis  of  the  system  of  equations  resu  lts  in  the  eigenvalues.  A;,  and  eigenvectors,  The  generalized 
coordinates  are  then,  n  —  In  the  ease  of  vibrational  analysis,  the  resonant  frequency  of  the  fh  eigenvalue  becomes 
iDi  =  \/A/-  The  resulting  set  of  equations  are  normalized  by  the  mass,  such  that  <prM(p  —  L  The  normalized  mass, 
damping,  and  stiffness  matrices  are  now  diagonal  and  scaled  appropriately.  A  reduced  order  mode!  is  formulated  by 
truncating  the  mode  set  at  a  desired  level,  and  all  of  this  information  is  put  into  an  input  file  (Kestrel  .csd  file)  for  the 
CFD-CSD  simulation* 

For  the  purposes  of  this  paper,  the  simplest  structural  model  was  used;  the  wing  material  is  homogeneous,  thus 
the  mass  and  damping  matrices  are  the  identity.  The  wing  structure  is  modeled  by  a  finite  dimensional  rectangular 
beam  using  brick  elements*  A  finite-element  method  (FEM)  code  was  employed  in  MATLAB^  to  solve  for  the 
corresponding  eigenvalues  and  eigenvectors  of  the  linear  system.  Assuming  a  linearly  elastic  material  allows  for 
changes  to  the  material  properties  by  simply  changing  the  generalized  mass,  damping,  and  stiffness  which  are  part  of 
the  eigen  solution*  In  addition,  these  values  are  accessible  directly  in  the  Kestrel  structure  input  file,  which  even  allows 
for  changes  to  these  values  individually  for  each  mode.  However,  one  must  be  careful  to  not  exceed  the  deformation 
limitations  imposed  by  the  assumption  of  a  linear  model. 


C.  Open-Loop  Forcing 


To  augment  the  flow  features  over  the  surface  of  the  wring,  blowing  slots,  as  shown  in  Figs*  2  and  3,  were  integrated 
along  the  wing.  The  base  of  the  slot  was  set  as  a  source  boundary  condition  such  that  the  magnitude  of  the  mass  flow 
through  the  patch  could  be  adjusted  using  the  Kestrel  SDK  features  described  above*  The  actuation  strength  is  then 
governed  by  the  mass  flow  through  the  source  patch  at  the  base  of  the  forcing  slot.  To  normalize  this  mass  flow  with 
the  momentum  in  the  flow,  the  normalized  momentum  coefficient,  C^,  is  defined  as  follows, 


"  P^UlAref 


(3) 


Because  the  velocity  vector,  F,  is  set  normal  to  the  patch  this  value  becomes  the  magnitude  of  the  vector,  VjeJ . 
Also,  because  exit  velocities  are  strictly  less  than  M  -  0*1,  the  velocity  at  the  jet  orifice  is  incompressible;  thus,  the 
density  is  constant  throughout  the  flowfield  and  can  be  removed  from  equation  (3),  The  expression  tor  the  momentum 
coefficient  can  then  be  simplified  to 


r 

91  tew 


(4) 


l  he  area  ratio  of  the  jet  orifice  (Ajef  —  I0~5nr)  to  the  reference  area  {Ayef  =  0.06m2 )  is  equal  to  1 .67  x  10  4,  and  the 
velocity  ratio  for  these  preliminary  results  was  chosen  such  that,  Vjel  =  Q.5CL,  giving  a  momentum  coefficient  of  — 
4.17xi0  *  per  slot.  Since  4  slots  are  used,  the  total  momentum  being  induced  into  the  flow  field  is  1 .67  x  1 0  "4* 


III*  Results 

The  results  below  will  be  presented  in  stepwise  manner,  where  the  three  regimes  of  fluid-structure  interaction  as 
identified  in  Section  [  are  presented.  First,  results  are  shown  for  the  rigid  wing  for  both  unforced  and  open-loop  forced 
simulations.  This  provides  insight  to  the  control  efficacy  of  the  blowing  ports  on  the  rigid  wing.  Secondly,  the  structural 
model  of  the  NACA0G18  is  presented  which  includes  spatial  modes,  stiffness,  and  resonate  frequencies.  Finally,  the 
frill  CFD-CSD  simulations  are  presented  which  show  both  unactuated  flutter  in  comparison  with  the  actuated  flutter 
simulation.  The  efficacy  of  open-loop  control  of  the  fluid-structure  interaction  problem  is  investigated  using  a  single 
actuation  amplitude.  The  resulting  forces  and  deformation  of  the  unforced  simulation  are  shown  to  be  different  from 
the  forced  simulation*  While  not  yet  conclusive,  these  results  provide  a  convincing  proof  of  concept  at  flow  control 
can  he  used  to  alter  the  flow  field  in  a  way  that  changes  the  structural  deformation* 
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A.  Rigid  wing 

For  comparison,  the  solution  for  a  rigid  wing  will  be  presented  first.  Figure  4  shows  the  surface  pressure  in  color  and 
surface  streamlines.  As  expected  at  this  angle  of  attack,  the  flow  is  attached  and  no  significant  flow  structures  develop. 
In  Fig.  4b,  a  dose-up  of  the  flow  field  near  the  blowing  ports  is  shown.  The  results  indicate  that  there  is  a  negligible 
effect  of  the  blowing  ports  on  the  flow  structures  when  they  are  turned  off. 


PRESSURE 
*  69200 
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a)  Full  wing. 


PRESSURE 


Figure  4,  Flow  field  of  Liu:  unde  formed  wing,  no  forcing.  Pressure  in  Pa. 

Once  the  blowing  is  turned  on  (Cp  =  1.67  x  10  4),  vortices  develop  at  the  blowing  slots  and  travel  downstream, 
figure  shows  the  flow  field  for  the  rigid  wing  when  the  blowing  slots  are  activated.  The  structures  are  visualized 
using  the  Q=0. 1  iso-surface.  Note  that  in  Fig.  5,  the  iso-surface  is  cut-off  at  x/c=G.5  for  visualization  purposes.  Due 
to  the  development  of  these  structures,  the  surface  pressure  changes  as  shown  in  Fig.  fi.  Low  pressure  can  now  be 
observed  near  the  blowing  slot  exits,  and  the  footprint  of  the  structures  is  visible  downstream.  As  noted  earlier,  at 
this  angle  of  attack,  the  effect  of  blowing  on  the  rigid  wing  was  expected  to  be  minimal,  and  the  results  confirm  this 
assertion.  The  main  purpose  of  these  simulations  was  to  provide  a  baseline  for  comparison  with  the  flexible  wing. 


PRESSURE  PRESSURE 


a)  Full  wing.  b)  Zoomed  view  near  blowing  ports. 

Figure  5.  Flow  field  of  (he  imddbrrnr  d  wing,  with  furring,  hosurfaee  of  normalized  Q=0.1.  Pressure  in  Fa. 

To  compare  the  effect  of  blowing  on  this  rigid  wing,  table  1  summarizes  the  lift,  drag,  pitching  moment,  and  roll 
moment  coefficients  of  the  wing  without  and  with  blowing.  As  shown  in  the  table,  blowing  increases  the  lift  coefficient 
by  approximately  7  percent,  which  in  turn  increases  the  magnitude  of  the  pitching  moment  and  the  rolling  moment 
due  to  the  changes  in  lift  occurring  in  the  outboard  region  of  the  wing  near  the  leading  edge.  No  measurable  effect  on 
the  drag  coefficient  was  noted. 
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PRESSURE  PRESSURE 


a)  Full  wing.  b)  Zoomed  view  near  blowing  ports. 


Figure  6*  How  Held  of  the  imrirformed  w  ing,  with  forcing*  Pressure  in  P«, 


fable  h  Coefficients  for  the  rigid  wing  with  and  without  blowing. 


cl 

CD 

c.m 
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0.430 

0.025 

0.025 

-0.075 

-0.085 

- 1 ,08 
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B.  Structural  modal  deformation 

As  outlined  above,  a  linear  modal  model  of  the  structure  was  developed.  The  cantilevered  NACA001S  wing  was 
assumed  to  have  a  rectangular  cross  section  with  matching  areas  for  simplicity.  For  the  finite  element  analysis,  beam 
elements  were  used  with  a  total  of  50  elements  discretizing  both  the  chordwise  and  the  span  wise  directions.  The 
thickness  was  discretized  by  a  single  element.  At  this  spatial  resolution,  convergence  was  achieved  and  the  results 
matched  experimental  static  deformation  measurements  and  the  analytical  solutions  of  the  modal  frequencies*  For 
the  coupled  fluid-structure  interaction  simulations,  a  5-mode  model  was  chosen  which  included  the  first  and  second 
spanwise  bending  modes  (Fig*  and  b),  an  in-plane  bending  mode  (Fig.  7c),  and  the  first  torsional  mode  (Fig.  7d). 
The  final  mode  in  the  model,  which  contributed  negligibly  to  the  deformation,  was  the  third  bending  mode. 


S  ^ 

a)  First  bending  mode  b)  Second  bending  mode  c)  In-phtne  bending  mode  d)  First  torsional  mode 

Figure  7*  The  first  four  modes  of  deformation  of  the  rectangular  beam. 


The  structural  parameters  of  the  wing  were  tailored  based  on  modal  parameters  to  approach  a  fi otter  limit  cycle 
oscillation*  The  damping  and  torsional  rigidity  was  estimated  through  a  flutter  prediction  technique  and  companion 
experiments  to  develop  a  wing  structure  which  naturally  oscillated  at  the  given  Mach  number  and  dynamic  pressure. 
[  J  A  summary  ot  the  structural  parameters  is  gi%ren  in  Table  7*  Note  that  the  bending  stiffness  is  much  larger  than  the 
torsional  stiffness  to  excite  the  torsional  instability  (he*,  a  limit  cycle  oscillation  in  angle  of  twist)  commonly  known  as 
stall  flutter  as  opposed  to  classical  flutter  in  which  the  torsional  and  bending  modes  couple  to  produce  a  pitch-plunge 
limit  cycle  oscillation  (LCO)* 

With  this  structural  model,  the  simulations  described  above  for  the  rigid  w  ing  were  repeated  for  the  elastic  w  ing  to 
assess  the  changes  in  the  forcing  effectiveness  due  to  the  deformation  of  the  wing. 
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labte  2.  Structural  parameters  for  the  aeroelastic  s  Limitations 


Bending: 

K„  [N/m] 

;  [N/m-s] 

m  [kg] 

K 

freq  [hz] 

64412 

2.84E-02 

0.S  1 

460085 

107,9 

Torsion: 

Ke  [N-m/rad] 
2.00 

n  rkg-m-s2/radl 
2.84  E-02 

i  [kg“in2/rad] 
2,51  E-04 

K 

7963 

freq  [hz] 
14.2 

Firsts  the  lime  history  of  the  lift  and  drag  coefficients  are  plotted  in  Fig.  and  the  time  history  of  the  wing  tip  twist 
and  bending  displacement  is  shown  in  Fig*  9.  Results  for  both  the  unforced  (C^  —  0}  and  forced  (C^  =  1 .67  x  10  4) 
aeroelastic  simulations  are  plotted*  As  shown  in  Fig*  8a,  the  oscillations  in  the  lift  coefficient  are  very  effectively 
reduced  by  the  blowing  because  it  induces  separation  at  lower  angles  of  attack  than  in  the  unforced  case*  Interestingly, 
after  the  initial  transient,  as  the  tip  twist  reduces  due  to  elastic  forces  (t  >  0.06st  see  Fig*  9a),  the  results  with  blowing 
show  a  large  reduction  in  the  amplitude  of  the  twist  oscillation  even  though  the  mean  lift  coefficients  does  not  seem  to 
be  reduced  significantly.  As  observed  for  the  lift  coefficient,  the  drag  coefficient  obtained  from  the  forced  simulation 
also  shows  a  significant  reduction  in  oscillation  amplitude*  Finally,  in  Fig.  ’  ,  the  wing  tip  deflection  due  to  bending 
shows  significant  changes  as  well.  In  the  case  of  flow  control*  the  initial  deflection  is  slightly  reduced,  but  as  the  twist 
angle  reduces,  the  bending  deformation  oscillates  around  its  new  mean  at  a  significantly  reduced  amplitude  and  at  a 
frequency  that  is  approximately  four  times  larger  than  in  the  unforced  case* 


Time  fs]  Time  [s] 


a)  Lift  cutflkknl. 


b)  Drag  coefficient. 


Figure  8*  Force  coefficients  for  unforced  {red)  and  forced  (blue)  simulations. 


To  better  illustrate  the  effectiveness  of  blowing  on  the  aeroelastic  behavior  of  the  flexible  wing,  a  phase  plot  is 
shown  in  Fig.  In,  plotting  the  lift  coefficient  against  the  tip  twist  angle  of  the  wing.  The  unactuated  case  (red)  shows 
an  initial  transient  due  to  the  change  from  rigid  wing  to  flexible  wing.  This  transient  is  due  to  the  onset  of  twist 
deformation  and  thus  generates  more  lift.  At  large  enough  angle,  the  flow  separates  and  the  elastic  forces  in  the 
structure  reduce  the  angle  of  twist.  For  the  unforced  simulation,  a  slightly  damped  LCD  results  and  an  orbit  about  an 
equilibrium  point  is  observed.  In  contrast,  the  actuated  simulation  (shown  in  blue  in  Fig.  10)  displays  the  same  initial 
transient  with  a  slightly  reduced  amplitude  (see  also  Fig,  9a)*  More  importantly,  the  phase  plot  clearly  shows  that 
the  actuated  wing  reaches  the  equilibrium  point  much  quicker  than  the  unactuated  wing*  From  a  dynamical  system 
perspective,  the  actuation  on  the  flexible  wing  has  increased  the  overall  damping  in  the  aeroelastic  behavior. 

Figure  1 1  shows  a  snapshot  of  the  instantaneous  fiowr  field  at  the  end  of  each  run,  t=Q.3s.  As  shown  in  the  figure, 
even  in  the  unforced  case,  Fig.  1 1  a,  the  control  ports  have  an  effect  on  the  flow  field  in  the  outboard  section  of  the 
wing*  The  data  indicates  that  the  flow  field  separates  further  upstream  behind  the  control  ports  when  compared  to  the 
inboard  section  of  the  wing.  In  contrast,  when  blowing  is  turned  on.  Fig*  !  lb,  the  skin  friction  coefficient  increases 
significantly  locally  around  the  blowing  ports,  indicating  the  development  of  streamwise  vortices*  which  result  in 
attached  flow  and  a  slightly  higher  drag,  which  may  be  partially  offset  by  the  reduced  skin  friction  inboard  of  the 
blowing  ports. 

When  forcing  is  turned  on,  the  flow  structures  identified  by  the  Q-criterion  (Fig*  13)  change  significantly  from 
the  unforced  case*  The  result  of  these  structures  generated  by  the  blowing  input  is  a  reduction  in  lift,  increase  in 
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Figure  9,  Tip  deformation  fur  im forced  (red)  ami  forced  (blue)  simulations. 


Figure  10,  Phase  plot  of  the  lift  coefficient  verse  the  tip  twist  for  the  unforced  (green)  and  forced  (blue)  simulation. 


a)  Unforced  h)  Forced 

Figure  11.  Instantaneous  deformed  geometry  and  skin  friction  coefficient  for  the  unforced  simulation.  t-fOs 
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drag,  and  a  reduction  in  the  magnitude  of  the  rolling  moment  (see  table  3).  This  is  almost  a  complete  reversal  of  the 
trends  seen  for  the  rigid  wing,  where  blowing  increased  the  lift  and  the  magnitude  of  the  pitching  moment  and  rolling 
moment.  Also  note,  comparing  tables  I  and  ??,  that  the  lift  and  drag  coefficients  arc  significantly  higher  for  the  flexible 
w  ing,  furthermore,  the  pitching  moment  is  reduced  and  the  rolling  moment  increased.  These  dramatic  changes  in  the 
performance  due  the  flexibility  of  the  wing  are  not  yet  fully  understood  and  are  currently  being  investigated. 

PRESSURE  PRESSURE 


a)  Full  wing.  b)  Zoomed  view  near  blowing  ports. 

Figure  T2*  Flow  field  of  I Tu?  deforming  wing  without  blowing. 


PRESSURE  PRESSURE 

.  W200  --  _  69200 


I  68000 
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a)  Full  wing. 


b)  Zoomed  view  near  blowing  ports* 


Figure  13*  Flow  field  of  I  he  deforming  wing  with  active  blow  ing. 


fable  3.  Coefficients  for  the  flexible  wing  with  and  without  blowing. 
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IV.  Conclusions 

In  this  paper,  HPCMP  CREAT  h  -AV  Kestrel  was  used  to  simulate  the  flow  around  a  rectangular  wing  with  open-loop 
flow  control  (steady  blowing).  Actuation  was  introduced  through  four  blowing  slots  in  the  outboard  section  of  the 
wing  near  the  leading  edge.  To  assess  the  efficacy  of  flow  control  to  alter  the  deformation  of  an  aero  elastically  soft 
wing,  the  flowfidd  around  a  rigid  wing  was  compared  to  the  one  around  the  deforming  wing.  The  results  indicate  that 
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control  by  blowing  has  a  pronounced  effect  on  the  flow  field  and  therefore  on  the  aerodynamic  coefficients.  The  lift  is 
increased,  as  are  the  pitching  and  rolling  moments. 

When  aeroelastic  deformation  is  considered,  the  picture  changes  significantly.  First,  the  flexible  wing  showed 
higher  lift  and  drag  due  to  the  increased  local  angle  of  attack  in  the  outboard  section  of  the  wing.  Furthermore,  the 
pitching  and  rolling  moments  were  significantly  reduced.  The  effect  of  these  changes  in  the  forces  and  moments  on 
the  overall  performance  of  the  aeroelastic  wing  have  not  been  completely  analyzed  because  the  time  history  acquired 
is  too  short.  Once  the  aeroelastic  wing  simulations  have  reached  their  dynamic  equilibrium,  further  analysis  of  the 
effect  of  forcing  will  be  performed. 
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